Robust control of individual nuclear spins in diamond 
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Isolated nuclear spins offer a promising building block for quantum information processing sys- 
tems, but their weak interactions often impede preparation, control, and detection. Hyperfine 
coupling to a proximal electronic spin can enhance each of these processes. Using the electronic spin 
of the nitrogen- vacancy center as an intermediary, we demonstrate robust initialization, single-qubit 
manipulation, and direct optical readout of ^^C, ^^N, and nuclear spins in diamond. These 
results pave the way for nitrogen nuclear spin based quantum information architectures in diamond. 
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The long coherence times of isolated nuclear spins have 
attracted considerable interest for applications in quan- 
tum information science [l|, 0, H, |3,H, H, 0, Hi , offering 
the potential for devices exhibiting coherent evolution 
over timescales exceeding several seconds Q. A signif- 
icant challenge in such applications is development of 
techniques to prepare, manipulate, detect, and couple 
individual nuclear spins. Some of the most promising ap- 
proaches use electronic spins as an intermediary to access 
single nuclear spins in solid state [1, 0, 0, S B IS, IHl 113 
or atomic [7, 8] systems. While atomic systems enable 
optical preparation and near-perfect isolation of nuclear 
spins, nanoscale solid-state systems offer the possibility 
of manipulating and coupling spins on fast timescales. 
Atom-like defects in solids ^ Q [H, d, [13, [Hi offer 
an intriguing interpolation between these approaches, al- 
lowing optical techniques in a setting amenable to fast 
control. 

We use the electronic spin of the negatively charged 
nitrogen- vacancy (NV), an atom- like defect, to access in- 
dividual nuclear spins in the diamond lattice. Even at 
room temperature, the NV spin is readily polarized and 
manipulated with optical and microwave fields [l^, yet 
it can exhibit millisecond coherence times [13, [lB| ■ More- 
over, its hyperfine interactions with nearby nuclear spins 
have allowed observation of single nuclear spin preces- 
sion [13, [T6l.[l8l. e ntanglement [Tsl . [Tsj , and two-bit quan- 
tum gates Hi As a result, the NV center has quickly 
become a strong candidate for spin-based quantum com- 
putation [Tsj and communication [20j applications. 

Three nuclear spin species commonly interact with the 
NV center, offering different properties that lend them- 
selves to different applications. Experimental manipula- 
tion of single nuclear spins associated with the NV cen- 
ter has emphasized / = \ ^^G [11, [H, [H, [11 [Ml and 
implanted I = \ ^^N nuclei [l^], while magnetic reso- 
nance in / = 1 ^''N has been observed in ensemble exper- 
iments [1^, [l^l- With the strongest hyperfine coupling, 
the nearest- neighbor ^■^C nuclear spins provide fast inter- 
actions but shorter coherence times [3, [3; like all ^^C 
nuclei, they occur randomly as isotopic impurities in the 
lattice. In contrast, the ^^N and ^^N species are weakly 
coupled, but they occur deterministically on the nitrogen 
of the NV center in natural (-'^'^N) and isotope implanted 



(^^N) diamond [22,[2l]- Scalable schemes - based, for ex- 
ample, on arrays of dipole-coupled electron-nuclear spin 
pairs or optically connected quantum registers [H [lol - 
are thus likely to favor the consistent properties of the 
nitrogen nuclear spins. Moreover, the longest electron 
spin coherence times have been observed only in highly 
isotopically purified, unimplanted diamond (itI . [l8 j, in 
which ^■^C and ^^N spins are exceedingly rare. These 
considerations point towards the development of control 
techniques that will function for ^^N and ^''N as well as 
^^C nuclear spins. 

In this paper, we demonstrate dynamic polarization, 
resonant manipulation, and direct measurement of sin- 
gle ^■^C, ^"^N, and ^^N nuclear spins. These techniques 
are robust in the sense that they do not rely on / = i 
spin properties or weak (selective) electron spin reso- 
nance (ESR), and give rise to signals as strong as the 
ESR signal. Moreover, these techniques share a common 
origin in the nonsecular components of the hyperfine in- 
teraction between the NV electron spin and the nuclear 
spin. The terms in the hyperfine interaction that give rise 
to electron-nuclear spin flip-flops allow robust polariza- 
tion , greatly enhance the Rabi frequency for nuclear 
magnetic resonance [llj, and enable direct optical read- 
out of the nuclear spin. These mechanisms apply to all 
species of nuclear spin coupled to the NV center, and we 
show that it may be possible to extend them to work in 
arbitrary magnetic fields. 

In our experiments, NV centers in CVD-grown high- 
purity diamond (Sumitomo, natural isotopic abundance) 
are isolated using confocal microscopy, allowing selection 
of NV centers with different proximal nuclear spins. The 
centers are illuminated by 4mW of 532nm light for polar- 
ization and fluorescence detection, while Helmholtz coils 
and/or a permanent magnet are used to apply magnetic 
fields between and 500G. Electron spin transitions are 
driven by resonant microwaves (MW) running through 
a 25 /im copper wire mounted on the sample; the same 
wire carries radio-frequency (RF) signals resonant with 
nuclear spin transitions. 

The NV center in diamond has an electronic spin S = 1 
in both its ground state (GS) and optically excited state 
(ES) [lil; both orbital configurations exhibit a zero-field 
splitting (A = 2.87 GHz (GS) or D = 1.42 GHz (ES)) 
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between the nis = and nis = ±1 spin projections along 
the NV axis (z) [1^. The electronic structure of the NV 
center also includes at least one long-lived spin singlet 
state of intermediate energy that preferentially depop- 
ulates the ES nis = ±1 states and decays into the GS 
rus = state. The singlet state provides the mechanism 
for polarization and detection of the electronic spin state: 
under optical illumination, the m^. = ±1 states fluoresce 
more weakly than the = states before being polar- 
ized into the to, = state 1191. 
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FIG. 1: (a) Ground state (GS) and excited state (ES) spin 
sublevels of the NV center, showing the excited state level 
anti-crossing (ESLAC) near 510G. Inset shows hyperfine lev- 
els and transitions for the ^''N isotope, (b-d) ESR for an NV 
center coupled to a proximal nuclear spin in B=65 G(top) 
and 509 G (bottom), (b) a nearby "G exhibits P = 91(5)% 
polarization; parentheses indicate the standard deviation sta- 
tistical error in the final digit, (c) ^''N: P = 95(3)%. (d) "N: 
P = 95(3)%. 

In certain magnetic fields, optical spin polarization can 
be extended to proximal nuclei [12, H^. Near 510G ||z, 
the ES TOs = — 1 and TOs = states become degenerate 
(see Fig. la), allowing electron- nuclear spin flip-flops to 
occur near-resonantly; repeated cycling through the sin- 
glet state thus results in polarization into = with 
maximal to/ (see Figlb-d). This mechanism was pre- 
viously used to polarize ^^7V and nearest-neighbor ^^C 
nuclear spins [22| and also polarizes the 1 — \ species 
^^N via the two-step process illustrated in Fig la. Non- 
nearest- neighbor ^•^C nuclei may also be polarized in this 
manner (see Fig. lb); however, not all "'^^C locations ex- 
hibit such polarization (data not shown) most likely be- 
cause their hyperfine axes are not aligned in the ground 
and excited states [l^, . 

Optical pumping at the excited state level-anticrossing 
(ESLAC) provides the polarization necessary to observe 
nuclear magnetic resonance (NMR) . Working within the 
TOs = — 1 manifold, we measure the nuclear spin projec- 
tion as a function of both RF frequency and duration 
(see Fig. 2), and fit our data to the expected form for 
square pulse excitation to extract the resonance frequen- 
cies ijJn and Rabi frequencies r2„ for all three nuclear spin 
species. Observed values of a;„ are consistent with pub- 
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FIG. 2: (a) Pulse sequence. Following polarization, a MW n 
pulse drives the system into the GS = — 1 state, where RF 
excitation excites nuclear spin transitions. To measure the re- 
sulting population in the nuclear spin levels, a second MW tt 
pulse flips the electron spin back into rris — 0, and subsequent 
fluorescence detection reveals the average nuclear spin projec- 
tion. All data is taken in a magnetic fleld of 509G \\z. Data 
for i^G, I'^N, and ^''N is shown in (b,e,h), (c,f,i) and (d,g,j) 
respectively, (b-d) Energy level diagrams for hyperflne levels; 
arrows indicate transitions driven by RF and MW; shading 
indicates the = — 1 manifold, (e-g) Fluorescence change 
((signal - reference) /reference) vs RF frequency; flt shows the 
expected form for a RF square pulse; the resonance frequency 
is shown with the statistical standard deviation for the flt. (h- 
j) Fluorescence vs RF pulse duration on resonance (frequency 
indicated on corresponding RF spectrum) superposed on a si- 
nusoidal flt. Short-time data (green) shows contrast of electron 
spin nutations for comparison. 



lished values [2J|, but the Rabi frequencies Qn greatly 
exceed the naive expectation J7„ = I'ynBnpl ^ a few 
kHz (we measure Brp ~ 5 — 10 G). Indeed, the nuclear 
spin nutation rates are dominated by a second order pro- 
cess akin to a chemical shift, involving an electron spin 
flip followed by an electron-nuclear spin flip-flop, with 
~ IuBrf + [HI. The enhancement oc- 

curs for all species of nuclear spin, enabling fast nuclear 
spin manipulation with low demands on RF power. 

The readout mechanism for the nuclear spin nutations 
shown in Fig. 2 is remarkably robust, yielding an NMR 
contrast as great or greater than the ESR signal - even 
for the nitrogen nuclear spins. Earlier nuclear spin read- 
out schemes have relied on nuclear-spin-selective MW 
pulse(s)within the GS to transfer nucle ar popu lations to 
measureable electron spin populations \l4 . fig , [TH . Such 
readout methods do not work as well for ^^N and ^*N 
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because their GS hyperfine interactions do not greatly 
exceed the electron spin dephasing rate. In contrast, by 
working at the ESLAC we can take advantage of the 
larger ES hyperfine parameters (A ~ 50 MHzj25l. [26j l 
to directly read out the nuclear spin state: Whereas 
the TUs = state with maximal nuclear spin projection 
fluoresces strongly, other nuclear spin projections pass 
through the dark singlet state during the polarization 
process, leading to reduced fluorescence. The reduction 
in fluorescence depends on the number of passes through 
the dark state required to reach the polarized state, so 
that the average nuclear spin projection can be directly 
measured with a contrast as great as or even exceed- 
ing the ESR signal. To illustrate the different fluores- 
cence levels. Figure 3b compares an ESR spectrum with 
TO/ = -1-1 (^''N) to one with to/ = (created using an 
NMR TT pulse on the transition indicated in Fig. 3a). 
The ESR spectra corroborate the nuclear spin projection 
and illustrate a significantly reduced fluorescence level 
for 771/ = as compared to mj = 1. While our mea- 
surements of this effect are aimed at direct nuclear spin 
readout, the same mechanism has recently been explored 
in great detail in (3]| for applications in enhanced elec- 
tron spin measurement. 
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FIG. 3: (a) Energy levels for GS mg = ^'^N nuclear spin 
manipulation, (b) Pulsed ESR spectra starting in ma = 0, 
mi = 1 (top) or ms = 0, mj = (bottom). Dashed lines 
indicate expected ESR positions for the three hyperfine lev- 
els, (c) Nuclear Rabi nutations within the ma = manifold 
(d) NMR spectrum obtained using the pulse sequence illus- 
trated below the data. A fit to the expected square-pulse 
lineshape yields resonant frequencies Wi — 5.094(1) MHz and 
u!2 = 4.7908(7) MHz. (e) Nuclear Ramsey fringe experiment. 
The downward trend is likely due to electron spin relaxation, 
which occurs at a rate Ti ~ 1 ms at room temperature [l9| . 

The preparation, control, and detection techniques de- 
scribed above enable straightforward measurement of nu- 
clear spin energy levels and coherence properties. As an 
example application, we measure some of the ^^N hy- 
perfine parameters to greater precision than current val- 
ues. Fig. 3d shows direct NMR signals for ^''N within 



the TOs = manifold. Combined with NMR resonant 
frequencies in the to^ = 1 manifold of w = 2.9373(2) 
and 6.9580(2) MHz (data not shown), we find A\\ = 
-2.162(2) MHz and P = -4.945(5) MHz, consistent with 
current published values [13, [3l[ . Coherence properties 
of individual nuclear spins are also readily obtained: for 
example. Fig. 3e shows a 2 kHz detuned Ramsey fringe 
experiment on the transition indicated in Fig. 3a. Al- 
though coherence properties for specific NV centers will 
vary with the mesoscopic distribution of ^'^C impurities, 
the long dephasing time T2 ~ few ms observed here bodes 
well for ^^N-based schemes. In short, these techniques 
open the door for complete characterization of individ- 
ual nuclear spins associated with the NV center. 

Resonant RF excitation can drive transitions in ar- 
bitrary nuclear spin species in any magnetic field, but 
the polarization and direct readout mechanisms function 
only in the vicinity of the ESLAC near 510G. This limita- 
tion is severe, especially for applications requiring optical 
interconnections between NV centers [l^l or stimulated 
emission-depletion techniques for sub-wavelength resolu- 
tion [131. Near the ESLAC, any light used to excite elec- 
tronic transitions can destroy information stored in nu- 
clear spins. While magnetic field switching techni ques 
developed in the context of fast field cycling NMR [23| 
may allow rapid transitions between two magnetic field 
regimes, careful examination of the NV ES level structure 
indicates an experimentally simpler alternative. 

Away from the ESLAC, the electron-nuclear spin tran- 
sitions that lead to optical nuclear spin polarization can 
be brought into resonance by microwave excitation within 
the ES manifold. However, under most conditions, the 
dominant effect of the microwaves will simply be to flip 
the electron spin, circumventing the hyperfine-induced 
flip-flops we wish to enhance. To induce the desired 
transitions, higher-order processes must prevail: For ex- 
ample, the axial (||z) component of the microwave mag- 
netic field and the nonsecular hyperfine terms together 
induce electron-nuclear spin fiip-flops through a second- 
order process, while transverse components (-L z) of the 
static and microwave fields yield similar results through 
a third-order process. 

To model microwave-induced polarization, we begin 
with the Hamiltonian for the NV-^^'A^ spin system in the 
excited state, 

H = DS^,+PI^, + (Bo + Bi cosujt)-{gfiBS - 7„I)-f AI-S, 

(1) 

where we have allowed for static and oscillatory magnetic 
fields and a quadrupole splitting P, and assumed isotropic 
g- factors and a contact hyperfine interaction A 26.]. We 
then use Floquet theory [201 to calculate the probability 
for different nuclear-spin-changing transitions in the ex- 
cited state as a function of the microwave frequency lu.Iii 
the case of ^''N, these probabilities can be incorporated 
into a simple 3-level rate equation model because sym- 
metry and the quadrupole moment conspire to keep the 
nuclear spin quantization axis \ \z for moderate magnetic 
fields. The transition probabilities are incorporated into 
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FIG. 4: (a) Theoretical positive {mi = mi = 1 plus 
mi = — 1 ^ mi = 0) and negative transition probabilities 
and equilibrium polarization for ^'*N in Bo = 505 G, p^bBi = 
20z MHz , with A = 50 MHz, P = 5 MHz, D = 1.42 GHz, 7„ 
= 0.3077 kHz/G, 7e = 2.799 MHz/G,7 = 2.5r, fee, = IQ-^F. 
(b) Same as for (a), but with 5a*sBi = 20x + 20z MHz 
and Bo = 40x + 48i G. (c) Polarization of in Bo = 
40a; + 48£ G. A 2/is pulse of MW and green light is followed 
by a 1/is delay before probing the nuclear spin populations 
in the ground state with weak ESR. Error bars indicate one 
standard deviation (statistical), (d) Polarization of ^^N under 
the same experimental conditions. 



enough to enable the third-order polarization process, we 
apply optical excitation and microwaves before taking a 
GS ESR spectrum. We fit the GS ESR signals with three 
(or two, for ^^N) Lorentzians constrained to have the 
same width, and use the depth of the fits to estimate 
the population in each of the nuclear spin projections. 
We thereby observe a small degree of polarization for 
both ^''N and ^^N nuclear spins (see Fig. 4c and d), and 
our results in Fig. 4c agree qualitatively with the model 
predictions in Fig 4b. 

Although our experiments are currently constrained by 
the orientation of Bi, our model indicates that precise 
orientation of Bq and Bi [3Ql] may allow full polarization 
of ^"'N (Fig 4a). Moreover, because each polarization step 
involves the metastable singlet state, such a polarization 
technique will also provide a direct nuclear spin readout 
mechanism. With careful engineering, it may be possible 
to enjoy the robust nuclear spin polarization and readout 
now available at the ESLAC at arbitrary magnetic fields. 

In conclusion, we have demonstrated initialization, ma- 
nipulation, and direct readout nuclear spins associated 
with the NV center in diamond, and proposed a mecha- 
nism to extend these techniques to other magnetic field 
regimes. Such control over nuclear spins has ramifica- 
tions for a broad range of applications in quantum infor- 
mation science. Moreover, these techniques enable mea- 
surements of nuclear spin parameters, interactions and 
coherence properties that may guide development of fu- 
ture quantum information processing devices. 



the model using excitation rates 7 (for transitions within 
the rus — manifold) and F (for transitions between 
the nis = and — ±1 manifolds). A depolarization 
rate k^q completes the model, which we use to predict 
the equilibrium nuclear spin polarization P+i — P_i [2^ . 
The details of this model will be published elsewhere. 
Figure 4a shows ^^N nuclear spin flip probabilities and 
equilibrium polarization for static and oscillatory mag- 
netic fields Bo and Bi oriented along z, while Figure 4b 
illustrates the same quantities for off-axis magnetic fields. 

The orientation of the microwave magnetic field Bi is 
constrained by the geometry of our experiments to be 
~45 degrees from the NV axis. Working in a sufficiently 
large axial static magnetic field Bq ^ to distinguish the ES 
TOs = ±1 levels, with a transverse component Bq.x large 
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